The spin excitations of the block-antiferromagnetic K .sFei 6 Se 2 
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We study the spin excitations of the newly discovered block-antiferromagnetic state in 
K0.8Fe1.6Se2 using an effective spin Hamiltonian suggested in the literature. Interestingly in ad- 
dition to the usual Goldstone mode, there exist three other "optical" spin wave branches. These 
spin excitations are the analog of "optical phonons" in crystals with more than one atom per unit 
cell. We compute the spin wave dispersion and dynamic spin structure factor, all of which can be 
measured by inelastic neutron scattering experiments. We also computed ordering moment reduced 
by spin wave fluctuations and the uniform spin susceptibilities. 

PACS numbers: 



The discovery of the K x Fe2- y Se2[I] system stirred up 
a new wave of excitement in the field of iron-based super- 
conductors. K 2 .Fe2_ l/ Se2 are isostructural to the "122"- 
family of iron pnictides, e.g. BaFe2As2- At ambient pres- 
sure they exhibit the highest superconducting transition 
temperature 33K among iron chalcogenidespQ. 

Very recently a neutron scattering experiment [5] 
showed that at the special composition x = 0.8, y = 0.4 
K0.8Fe1.6Se2 exhibits a novel structure and magnetic 
state. In particular for T < 578K the Fe vacancies 
(which cause the deficiency of the iron content from 2 
Fe per formula unit) form a y/E x \/5 superstructure in 



each Fe plane as shown in Fig. 1(a) The nearest Fe-Fe 
distance between the four irons in each y5 X v5 unit 
cell is slightly smaller than that between the neighboring 
cells [H [S]. Moreover for T < 559K the system becomes 
magnetic. The magnetic moments of the four irons in 
each V5 x y5 unit cell align ferromagnetically along the 
crystalline c-axis. These aligned "super-moments" stag- 
ger from cell to cell to form a block-checkerboard antifer- 
romagnetic pattern [2]. At low temperatures the ordered 
moment per iron is approximately equal to 3.3/is. Since 
the iron valence is found to be +2(1], and LDA calcu- 
lations find K0.8Fe1.6Se2 to be an insulator with an ap- 
proximate 0.6eV energy gap[31l3, this implies that each 
iron is in spin 2 state. 

In addition to the novel structural and magnetic or- 
dering, transport measurement has shown samples with 
composition close to K0.sFe1.6Se2 are superconductors 
with transition temperature in the neighborhood of 30K. 
Ref.[2j |4j advocates the microscopic coexistence of su- 
perconductivity and block-antiferromagnetism. Others 
raises doubt about this conclusion^. We will restrain 
from discussing superconductivity in this paper before 
the experimental situation is clear. In the following we 
shall focus on the magnetic property of K0.sFe1.6Se2- 

Given the novel magnetic order described above it is in- 
teresting to ask what is the nature of the spin excitations. 
In particular the block-antiferromagnetic order suggests 



in addition to the usual gapless Goldstone modes (which 
can exhibit a small gap in the presence of spin-orbit in- 
teraction), there should be "optical" spin wave branches 
analogous to the optical phonons in crystals with more 
than one atom per unit cell. Given the fact that the or- 
dering moment is relatively large, the linearized Holstcin- 
Primikov Hamiltonian should be a very good approxi- 
mate description of the low energy spin dynamics. 

In Ref.[3. an effective spin Hamiltonian has been used 
to fit the LDA results. The Hamiltonian involves intra- 
block nearest and second neighbor interactions J± , J2 and 
the inter-block nearest and second neighbor interactions 
Jy , J2' • Because the iron magnetic moments are ob- 
served to align along the c-axis (perpendicular to the iron 
planes) [2j , there should presumably be a spin anisotropy 
term in the Hamiltonian. The minimum Hamiltonian we 
consider reads: 
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FIG. 1: (Color online) (a) The schematic representation of 
the vacancy ordered Vo x v5 lattice structure in each iron 
plane. Circles represent irons and the solid and dashed lines 
represent different types of spin interaction, (b) The ordered 
spin configuration according to Ref . [2] . The spins point along 
the c-axis (out of plane). 
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where A is the spin anisotropy, and 



-JvT\ 



J 2 r 



j 2 ,t% with i* (6 



1, T, 2, 2') being the adja- 



cency matrix defined by 



pt, _ f 1 : form the shortest 6-type bounds, 
lJ 1 : otherwise. 

(2) 

The bond types are illustrated in Fig. |l(a) The classi- 
cal ground state has the block-AFM order as shown in 



Fig. 1(b) For later convenience, the sign factor Q = ±1 
for i G the up(down)-spin block is defined here. 

In the following we assume each spin has S = 2. To 
study the spin wave excitations we use the linearized 
Holstein-Primakoff (HP) transformations Sf — Q(S — 
a\a,i) and Sf = y/2S(a\ +a,)/2. Using the magnetic unit 
cell with 8 Fe atoms, the spin Hamiltonian can be written 
in terms of HP bosons in the momentum space as 



H =2 S 



(3) 



where <\>\ a = {a\ a ,a 



qa) and 

Jibuti) + J2T 2 a p(q) + J o 8 a0 
<j\. 



(TO 



(4) 



Here a, (3 label the 8 sites in the unit cell, Jo = 2A — 2 J\ + 
Jy — J2 + 2./2', and <7i denote the Pauli matrices (do being 
the 2x2 identity matrix) . The Fourier transform of the 
adjacency matrix reads 



(g) 



rb iq-(ri-r 3 ) 



(5) 



Here denotes the position vector of site i, and q € the 
magnetic Brillouin zone (MBZ). Use the the Fe-Fe bond 
directions as x and y and set nearest Fe-Fe bond length to 
unity, the MBZ is spanned by the reciprocal lattice vec- 
tors Qi = (3tt/5, tt/5) and Q 2 = (-tt/5, 3tt/5). The high 




FIG. 2: (Color online) (a) The magnetic Brillouin zone, (b) 
The spin wave spectrum. The black lines represent the dis- 
persion in the absence of the spin anisotropy. The red lines 
represent the dispersion in the presence of spin anisotropy. 
As expected, spin anisotropy gaps the Goldstone mode. 



symmetry points are defined as T = (0,0), X' = Qi/2, 
Y' = Q2/2, M' = (Qi + Q 2 )/2 as shown in Fig.|2(a) 



For each q we perform the Bogoliubov transformation 
on the 16 x 16 matrix SJ(q) to obtain the spin wave op- 
erators and dispersions. This can be achieved by solving 
the generalized eigenvalue problem 







SJ a p{q) 



Vnq/3 



nq&Z 



(6) 



where Q nq > are the spin wave dispersions and 11 = 
1, • • • ,8 denotes the band index (keeping only the pos- 
itive eigenvalue solutions). Note that 03 is needed in 
Eq.(|6]) for Bogoliubov transformation of bosons. The 



wave functions are normalized to ^2 a (\u r , 



«| 2 ) 



1. The corresponding spin wave creation operator is 
given by b\ q = Y, a u nqa a \ a + v nqa a-qa- Due to the 
Goldstone theorem and the residual U(l) rotational sym- 
metry of the spins, we anticipate two branches of degen- 
erate gapless spin wave modes (which would be slightly 
gapped by the spin anisotropy). However, since there 
are eight spins per magnetic unit cell, there must remain 
6 gapful collective branches in the spin wave dispersion. 
These gapped spin wave modes are like optical phonons in 
the crystal with multi-atom per unit cell. These gapped 
spin waves are, in general, twofold degenerate due to the 
residual U(l) symmetry. 

To explicitly demonstrate the above anticipation in the 
absence of the uniaxial anisotropy (A = 0) we adapt 
the parameters J\ — — 4meV, Jy — lmeV, J2 = 3meV, 
J2' = 12mcV. Aside from J 2 the above values are given 
in Ref.(3] by fitting to LDA calculations. The reason 
we slightly reduced J2 is because using the value given 
in Ref.[3] we find the spin wave excitations are unsta- 
ble. For the chosen parameters, the resulting spin-wave 
dispersion is shown as the black lines in Fig. |2(b)| As ex- 
pected, there are two degenerate Goldstone modes. The 
rest of the spin wave excitations form three doubly degen- 
erate branches. In the following we shall refer to them as 
the "optical" spin waves. To demonstrate the effect of the 
uniaxial anisotropy we set A = 0.08meV and illustrate 
the corresponding spin wave dispersion using the red lines 
in Fig. |2(b)| The gap of the Goldstone mode at q = is 
found to be 2S' V /A(A + J v +2J 2 ') ~ 6meV. The optical 
modes are not much affected by the anisotropy as long 
as A is much smaller than the magnitude of all coupling 
constants. 

In table [EJ we list the analytical expressions for the 
optical spin wave excitation energies at T and M' (for 
A = 0). (Analytical expressions at other momentum 
points are not known to us.) In conjunction with the 
neutron data they can be used to determine the effective 
coupling constants J\ t 2 and Jx',2'- 
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TABLE I: Table of spin wave excitation energies (A = 0) 



Momentum 




Degeneracy 


r 




2S y /{J 1 + J2- J 2 >)(Ji - Jv +J2- J a >) 
2 V / 2S' v /(2Ji - Jv) (Ji - Jv) 


2 
4 
2 


M' 


2S^/Jl - (J v - 2J 2 + 2J 2 ,)J 1 + J 2 (J 2 - 2J 2 ,) + J V {J 2 , - J a ) 
25| Ji + ^/(Ji - J 2 ,)(Ji - Ji/ - Jv)\ 
2S\ Ji - v/(Ji - Jy)(Ji - Ji/ - J 2 ')l 


4 
2 
2 



Due to quantum and thermal fluctuations, the mag- 
netization \(Sf)\ is reduced from S to \(Sf)\ = S - SS, 
where the correction SS is given by 



SS 



d 2 q 



~,{ a qa a qot) l 



qGMBZ 

+ |w„q Q | 2 (n B (^nq) + 1) 



where ub is the Bose-Einstein distribution function and 
there is no summation over a. The moment correction 
SS is found to be the same on all sites. For the parame- 
ter Ji = -4, J v = 1, J 2 = 3, J 2 < = 12, A = 0.08meV, we 
find SS = 0.32 at zero temperature. This amounts to a 
magnetic moment per iron equal to 2 x (S — SS) = 3.36 
Bohr magneton (ub), in good agreement with the exper- 
imental deduced value. The temperature dependence of 



the ordering moment (per iron) is shown in Fig. 4(a) 



To the same leading order approximation, the zero- 
temperature dynamic structure factor is given by 



S±(q,u>) 



, A nq S(u - Q nq ), 



(8) 




T M' 2M' 3M' 4M' 5 M' 

FIG. 3: Dynamic spectral function along the cut through 
the F-M' line. pM' denotes the momentum point of 
(jJ7r/5, 2piv/5). 2M' and 4M' are equivalent to F point in 
the MBZ, while 3M' and 5M' are equivalent to M' . The spin 
wave dispersion is plotted using the spectral weight of the 
dynamic spectral function as darkness. Darker line implies 
greater spectral weight. 



where the spectral weight A nq reads 
A nq = ^Kqa + v* nqa ){u nq0 + v nq p)e lq < r «- r ?\ (9) 

Note that A nq is not periodic with reciprocal magnetic 
unit vector Qi and Q2 even though spin wave energy 
il nq is. This is similar to the matrix element effect in 
angle-resolved photoemission experiments. Fig. [3] shows 
the dynamic structure factor along the T-M' cut. In in- 
elastic neutron scattering experiments, the T-point "opti- 
cal" modes will be observable at (2n/5, 47r/5) point, and 
the M'-point "optical" modes will be visible at (tt, 0) 
point. 

It is also straightforward to compute the uniform sus- 
ceptibility. The transverse uniform susceptibility reads 



Xx = ^ Hm Y 
A 128 g^o ^ n 

n,a,p 



n 



(10) 



which is independent of temperature in the low temper- 
ature regime. Using our parameters, it can be evaluated 
to xx = 2.49/i^eV" 1 . The longitudinal uniform suscep- 
tibility is given by 



XII 



Ejl \ " 

256 ^ 

n,q L Ck 



y(a(\u nqa \ 2 + \v nqa \ 2 ) 



n' B (n nq )in) 



where where £ a = ±1 for a residing in spin up or down 
block respectively, and n' B denotes the derivative of Bose- 



fJ-/t*B 
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FIG. 4: Temperature dependence of (a) the per-iron magnetic 
moment, and (b) the longitudinal uniform susceptibility. 
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Einstein distribution function. Fig. |4(b)| shows the tem- 
perature dependence of x\\ ■ Below the temperature scale 
of the spin gap (~50K), x\\ exhibit exponential temper- 
ature dependence. While above that temperature scale, 
the linear-T behavior x\\ °^ T can be observed. 

In conclusion we have studied the spin wave excita- 
tions of the block-antiferromagnetic K0.sFe1.6Se2 using 
an effective spin model suggested by LDA calculation. 
There are both the Goldstone modes and the "optical" 
spin wave branches. After the spin wave correction, the 
ordering moment is found to be « 3.4/ljb per iron at low 
temperatures, which is close to the observed value. The 
observation of the optical spin wave modes will constitute 
a spectacular supporting evidence for the block antifer- 
romagnetism. 

Currently we do not have a good understanding of the 
block-antiferromagnetism and its relation to the struc- 
tural distortion from microscopic point of view. Under- 



standing it can shed light on the local electronic structure 
and correlation for the newly discovered K.;;Fe2-ySe2 and 
other pnictides. 
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